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We study the nonlocal spin transport in a lateral spin valve with multiple ferromagnetic (FM) 
electrodes. When two current-injecting and two spin current-detecting electrodes are all ferromag- 
netic, the number of possible nonlocal spin signal states is four at maximum. In reality, this number 
is reduced, depending on the inter-probe distance and the relative magnitudes of the spin resistances. 
Our theoretical results are in agreement with recent experiments of spin injection into an Al island, 
a carbon nanotube, and graphene. 
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I. INTRODUCTION 

Spin injection 1 ^ from a ferromagnetic (FM) metal into 
nonmagnetic (NM) materials (metal, semiconductor, in- 
sulator) is very important for device applications and for 
academic interest. The operation of the spin valve, a hy- 
brid structure of FM metal/NM/FM metal, depends on 
efficient spin injection from one FM electrode into the 
NM layer and spin detection in the other FM electrode. 
Typical examples of a spin valve include giant magne- 
toresistance (GMR) devices^, magnetic tunnel junctions 
(MTJ) 4 A&, FM/NM/FM nanopillars 7 ., etc. The mag- 
netoresistance has been enhanced up to a few hundred 
percent in the recent MgO-based MTJs^ and GMR and 
MTJ devices are already in commercial markets for mag- 
netic read heads. MgO-based MTJs are now being used 
for magnetic random access memory devices. However, 
these vertical spin valves have some difficulties in inte- 
grating them into semiconductor electronics. 

In order to integrate spintronic devices into semicon- 
ductor electronics, lateral spin valves are more desirable 
for, e.g., multi terminal devices. Due to the increased 
distance between terminals in lateral spin valves, the 
spin signal is more suppressed compared to vertical spin 
valves. Furthermore, spin injection and detection experi- 
ments in the two terminal geometry of lateral spin valves 
are complicated and obscured by other effects like the 
anisotropic MR, the anomalous Hall effect, etc. In order 
to avoid these undesirable effects, a nonlocal spin valve 
geometry 1 in which the spin current path is spatially sep- 
arated from the charge current path, was adopted, and 
spin injection and detection were clearly demonstrated 
in Al wiresi^. The nonlocal spin injection technique was 
also used to observ o 11 : 12 the (inverse) spin Hall effect in 
diffusive nonmagnetic metals (e.g., Pt with a spin diffu- 
sion length of 10 nm), which are characterized by a rather 
strong spin-orbit interaction. In these experiments, the 
separation of charge and spin currents, as well as efficient 
spin injection, is essential to observing the Hall voltage 



induced by the spin current. 

If significant magnetoresistance (MR) or nonlocal spin 
signals are to be achieved in lateral spin valves, the spin 
polarized current should be able to pass through the in- 
tervening nonmagnetic layer without losing or degrading 
too much of its spin polarization in the NM layer or at 
the FM-NM interface. Nonmagnetic materials with long 
spin diffusion lengths (SDLs) are most desirable and are 
required for successful operation of lateral spin valves. 
In this respect, recent experiments of spin injection into 
carbon systems, carbon nanotube^ and graphen o 14 ' 15 , 
are very intriguing. A large spin polarization and a large 
SDL are observed in carbon systems. In these experi- 
ments, the two voltage probes and thr two current probes 
were all ferromagnetic in contrast to previous experi- 
ments. Motivated by these experiments, we study in this 
work the spin transport in a nonlocal spin valve geome- 
try with four FM electrodes, as shown schematically in 
Fig. [H based on the one-dimensional spin drift-diffusion 
(SDD) equations. 



II. FORMALISM: SPIN DRIFT-DIFFUSION 
EQUATION 

Spin transport in spin valves can be understood 
theoretically based on the spin drift-diffusion (SDD) 
equations 1 ^, which are a reduced version of the spin- 
dependent Boltzmann equation. The SDD equations 
were shown to be valid when the mean free path is much 
less thaniS or comparable toU the spin diffusion length. 
The SDL is the length scale over which electrons can 
preserve their spin information. The finite SDL in the 
samples is caused by the spin-flip scattering due to the 
spin-orbit interaction, magnetic impurities, etc. Though 
the SDD equations are phcnomcnological, they have been 
very successful in explaining the main features of exper- 
iments qualitatively. 

The SDD equations are written down for the spin- 
dependent electrochemical potential fi a and an electric 
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current density j a . Here, a = ± represents the spin-up 
(+) and the spin-down (— ) states, respectively. In the 
one-dimensional device structure, the SDD equations can 
be written down in a matrix form as 



dx 2 



Ja 




(1) 
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Here ,D a is the diffusion constant for spin direction a — 

±, and T_| is the average spin-flip time for an electron 

from the spin direction + to — . a a is the conductivity 
for electrons with spin a, and e is the absolute value of 
the electron charge. Analyzing the eigenvalues and the 
eigenvectors of the matrix in the SDD equation, we can 
find the general solution of the SDD equations for the 
electrochemical potential and the corresponding current 
density^: 



As we shall show below, the use of properly defined ma- 
terial parameters highly simplifies the algebra. 

We consider the model spin valve structure in Fig. Q] 
The spin-polarized current I flows from F2 into N and Fl. 
The spin-dependent electrochemical potential (fJ.i±) and 
current (Ii±) in each FM electrode can be parameterized 
as 



-y S. 



-y S. 



Ui 



<j\A x a 2 A 2 
-v/x, ( + 



(6) 



(7) 



h 






X 








.a 


1 






1 


j 




a 



A 



BX 



± E. e -NA 

(7-1 2 



-\x\/X 



-1 



-i (3) 



(4) 



Here, a = a+ +c_ is the total conductivity, and A and B 
are parameters to be determined by the boundary con- 
ditions. The total or charge current density is constant 
and uniform in space: j = j+ + The spin diffusion 
length A is defined by the expression 
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In a one-dimensional device structure, it is more conve- 
nient in algebra to use the current instead of its density. 

(a) Fl 



The index i runs from 1 through 4, labeling the four FM 
electrodes. As will be clear below, Ui = Vi + A<Jj ; i and 
the shift A is introduced to take into account the current 
flow in the common electrode. That is, the electrochem- 
ical potential in Fl is shifted up with respect to that in 
F2, F3, and F4. Vi is the voltage drops in each FM elec- 
trode Fi, which is induced by the spin accumulation and 
diffusion. Ii measures the spin current leaking into Fz 
and should be determined by the Kirchoff rules at the 
junctions. The spin-dependent current is determined by 
the equation 
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Note that the i-th FM lead is contacted to the common 
electrode at x = Ai, A,, <7j, and /3j are the cross- 
sectional area, the spin diffusion length, the conductivity, 
and the bulk spin polarization in the conductivity of the 
i-th FM lead, respectively. IZi is the resistance of the FM 
lead over the spin diffusion length and is defined by the 
relation 
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FIG. 1: Schematic display of the lateral spin valve system 
with four ferromagnetic electrodes, (a) Ferromagnetic elec- 
trodes are labeled as Fl, F2, F3, and F4 from left to right. 
The base electrode is denoted as N. The current I is injected 
from F2 into N and drained to Fl. The voltage induced by 
the spin current is measured between F3 and F4 electrodes, 
(b) The one- dimensional model geometry of the spin valve 
system in (a). 



The spin-dependent conductivity <7j± can be written in 
terms of the spin polarization /3, as 



&i± = ^i{l±l3i). 



(10) 



In the common base electrode, the spin-dependent elec- 
trochemical potential ) and current {I± ) can be writ- 
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ten as 
e \M- 
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Note that the electrochemical potential shift, A, is de- 
fined by the relation 



(13) 



Due to the current flow between Fl and F2, the elec- 
trochemical potential is shifted up in the region x < L x 
with respect to the region x > L 2 . The spin-dependent 
current is determined by the equation 
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A, A, and a are the cross-sectional ell Cel. the spin diffu- 
sion length, and the conductivity of the base electrode 
N, respectively. 1Z = X/aA is the resistance of the base 
electrode, which is defined over the spin diffusion length 
of the base electrode. 

There are three sets of unknown parameters, V*, L, and 
Ji's, which we are going to determine by using the bound- 
ary conditions or the Kirchoff rules at the junctions. The 
above expressions of the electrochemical potential are 
constructed such that the charge current at each junction 
is conserved and flows from the F2 lead into the com- 
mon base electrode and finally into the Fl lead. Though 
the spin- flip scattering in the bulk is taken into account, 
any possible spin-flip scattering at the interface between 
the FM leads and the nonmagnetic base electrode is ne- 
glected. In this case, the spin current (I s =/+—/_) is 
also conserved at each junction and leads to the following 
constraint: 
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The electrochemical potential for each spin direction 
should satisfy Ohm's law at the junctions. At the i- 
th junction, the drop in the electrochemical potential is 
given by the expression 
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Ohm's law at the junctions gives the following relation: 

n ti ±i i± , (17) 
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where lZu± is the spin-dependent junction resistance at 
the interface between the base electrode and the z-th FM 
electrode Fi, and can be defined in terms of the spin 
polarization 7$ of the junction resistance as 



2-Ktj 
l± 7i 



(18) 



Ii± is the spin-dependent current passing through the 
i-th junction: 

h± = ^(l±p 1 )6 iil - I -(l±fa)8 ii2 ± I j. (19) 

In our work, two types of spin polarization are intro- 
duced. One is /3j, measuring the spin polarization in the 
bulk conductivity in the FM electrode Fi. The other one 
is 7i, which measures the spin polarization in the spin- 
dependent junction resistance between N and Fi. In our 
convention, (3 and 7 are positive when spin-up electrons 
are in the majority band while they are negative when 
the magnetization orientation is reversed or the spin-up 
electrons are in the minority band. 

For the algebraic manipulation, it is more convenient 
to introduce new parameters for the resistance: 
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and R = 1Z. With these new notations, we have 

V t = ±(1T 0i)IiRiTRYl A ij J j 

3 

+Rti(l T 7*) [/(I ± - 1(1 ± P 2 )h, 2 ± #1) 

Adding and subtracting the two equations, we find 

Vi = -('nR t i + 0iRi)ii + (i-0i'n)Rtii Si,i 

-(1-/3272)^2/^,2, (22) 

and 

(Ri + Ru)Ii — AijJj 
3 

= -(fo - j!)R n I + (f3 2 - 72)^2/ S it2 (23) 

For the algebraic manipulation, it is much more conve- 
nient to use the vector and matrix notations to rewrite 
the equations obtained from the boundary conditions as 

\V> = -(7Rt + /3R)|/>+(l-/?i7i)#ti^ |1 > 

-(1-/3272)^2/ |2>, (24) 
|/> = -2\J>-0 1 I\1>+0 2 I\2>, (25) 



and 



(R + Rt)|J > -RA\J > 

= -(/?! - 7 i)i?ti/|l > +C82 - 72)^ t 2/|2 > (26) 
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R and Rt are diagonal matrices with their diagonal el- 
ements Ru = Ri and R tii — R ti , respectively. We used 
the vector notations, for example, < I\ — (ii, I 2 , 13, h) 
and < 1| = (1,0,0,0). After some algebra, we find the 
expression of the desired electrochemical potential drops: 

^ = (-l) i+1 (R ti +^R i )5 i < 2 
2 

where the matrix G is defined by the relation 

-1 -1 



R + R, 



1 



RA 



(28) 



Note that the matrix G has the dimension of conductance 
and is independent of the magnetization configurations in 
the FM electrodes. The spin current Ii is also obtained 
as 



-Pita + flafo + GniPiRi + -yiRa) 
-G l2 {(3 2 R2 + i2Rt2)- 



(29) 



73 and Z4 measure the leakage spin currents into F3 and 
F4, respectively. Since Gij is independent of the mag- 
netization configurations, the signs and the magnitudes 
of the leaking spin currents do not depend on the mag- 
netization orientations of the F3 and the F4 electrodes 
(parallel or antiparallel relative to Fl and F2), but de- 
pend on those of Fl and F2. On the other hand, the 
voltage drops V3 and V4 change their sign when their 
magnetizations are reversed. 



III. RESULTS AND DISCUSSION 

In the previous section, we derived the voltage drops 
and leaking spin currents in the voltage probes F3 and 
F4. We can consider several different measurement ge- 
ometries from our analytic solutions. 

For the voltage probes F3 and F4, we have the spin 
leakage currents and the corresponding voltage drops 
from Eqs. (|2T|) and 



j- = G a (/3 1 R 1 +j 1 R t i)-G i2 ((3 2 R 2 +j 2 Rt2),(30) 



Note that the nonlocal voltage drop can be expressed 
as the product of the leakage spin current and the ap- 
propriately defined spin resistance. This spin resistance 
depends on the orientation of the magnetization. Our 
general results can be reduced to the simpler ones. When 
Fl is located very far away from F2 or when the distance 
between Fl and F2 or \L 2 — Ll\ is much longer than 
the SDL of the nonmagnetic electrode, our spin valve 
structure is reduced to that studied by the authors in 
Ref. [I?! where the effect of an additional FM electrode 
on the nonlocal spin signals was studied. Furthermore, if 
the distance between F3 and F4 is much longer than the 
spin diffusion length, our model spin valve is reduced to 
the conventional spin valve studied in Ref. [ToL 

The nonlocal spin signal in the model spin valve in 
Fig. [U is quantified by measuring the voltage difference 
between two FM electrodes, F3 and F4. Hence AV = 
V3 — V4 is the experimentally relevant quantity and the 
nonlocal spin signal is defined as R s — AV/I, where 

i=3,4 j=l, 2 

x( 7j R tj + fcRj). (32) 

This is the main result of our work. 

To get some insight about the nonlocal spin signal R s , 
let us start with an algebraic manipulation of the con- 
ductance matrix G. For the calculation of the inverse 
matrix G, we note the symmetry of G and G _1 can be 
rewritten as 



G- 1 



A 



D 



A T 
Tt B 

n \hR 

\flR T 2 

r 3 \hR 

\hR r 4 



Vi = -(jiRu + PiRi)Ii. 



(31) 



(33) 
(34) 
(35) 
(36) 

Here, ft = e -\L^+i-Li\/\ with i = i ? 2,3, and r l = R4 + 
Ru+\R with i = 1, 2, 3, 4. The 2x2 matrix T is real, and 
its hermitian T' is equivalent to its transpose T*. Though 
\T\ = 0, the determinants of A and B are nonzero so that 
we can find the conductance matrix G in a manageable 
form. The inverse of a 4 x 4 matrix is reduced to that of 
2x2 matrices, which is much simpler to compute: 



T = -f 2 R ( f l fl / 3 

r 2 \ 1 h 



_ ( A-^l-TB-^A- 1 ]- 1 -A^TB-^l-T^A^TB- 1 ]- 1 ^ 
<* - y-B- 1 T*A- 1 [l-TB- l T*A- 1 ]- 1 B^l - T^A^TB- 1 ]- 1 J' {67 > 

For the nonlocal spin signals such as ^3, 14 and V3, V4, only the off-diagonal block matrix is relevant. Since T contains 
the exponentially decaying factor with the scale of the spin diffusion length, we may approximate the matrix G in the 
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expression of a nonlocal spin signal as 



G31 G32 
G41 G42 



= -b-Wa-^i-tb-Wa- 1 ]- 1 « -.e^rU -1 

2(r 1 r a - i/fiil 2 )^ - i/fi? 2 ) U/afo - |i?)(r 3 - \R) f 3 (n - ±f?R){r a - J ' 



This approximate form is valid when the distance be- 
tween the current probe F2 and the voltage probe F3 is 
much larger than the SDL so that the matrix T is small. 
For a more accurate calculation, the exact form of the 
conductance matrix G is needed. However, the above 
approximate conductance matrix is enough to get some 
insight into the nonlocal spin signal R s . 

Since fi is an exponentially decaying factor over the 
spin diffusion length, the most dominant contribution in 
the nonlocal spin signal R s comes from G32. The nonlo- 
cal spin signal or the transresistance R s in Eq. (|3"2")) can 
be rewritten in decreasing order as 

R s = {isRts + &#3)G 32 (72-Rt2 + foRa) 
-(jsRts + /?3# 3 )G 3 i(7i^i + PiRi) 

-(j 4 Rt4 + PARA)G A 2{l2Rt2 + P2R2) 

+( 74 i?t4 + p4Ri)G u { 1 iR n + ftife). (39) 

Due to the exponential factors, the first line is dominant 
and the last line is the weakest. Obviously, the nonlocal 
spin signal will change its sign when the magnetization 
orientation of two neighboring FM electrodes F2 and F3 
is inverted. 

In experiments, the magnetic configurations are 
scanned by sweeping the external magnetic field along 
the geometrically parallel FM electrodes. With differ- 
ent widths of the FM electrodes, the coercive fields are 
different due to the different strengths of an easy axis 
anisotropy along the electrode direction. Starting from 
a large negative field, all the FM electrodes are aligned 
with the magnetic field. As the magnetic field is swept 
from negative to positive, the FM electrodes reverse their 
magnetization one by one when the magnetic field equals 
their coercive field. While sweeping the magnetic field 
from large negative to large positive, we have four dif- 
ferent magnetization configurations so that nonlocal spin 
signal will show four different states. When the magnetic 
field is now swept from large positive to large negative, we 
have again four different magnetization configurations. 
Under the time reversal operation or when the magneti- 



zation orientation of all four FM electrodes is reversed, 
we can readily see that the nonlocal spin signal R s does 
not change its sign. This means that there are only four 
different values of R s under the magnetic field sweeping. 
In principle, there can be 4 different states in the non- 
local spin signal, depending on the relative orientations 
of magnetizations. However, the number of actual states 
realized in experiments may well depend on the distance 
between the FM contacts, the relative magnitudes of the 
spin resistances in R s . We note that all four states were 
observed in Ref. H3 for a nonlocal spin valve with Al as 
a base electrode, three states in Ref. [H], and two states 
in Ref. [H 



IV. SUMMARY AND CONCLUSION 

In this paper, we studied the nonlocal spin transport in 
a lateral spin valve with four FM electrodes, out of which 
two are used as current injectors and the other two are 
used as spin current detectors or the voltage probes. We 
calculated the general expressions for the nonlocal sig- 
nals, such as the leakage spin current and the voltage 
drop due to the spin accumulation. Since there are four 
ferromagnetic electrodes in our model spin valve struc- 
ture, in principle, four different nonlocal spin signal states 
are possible. In real experiments, the number of observed 
states is found to depend on the inter-electrode distance, 
the relative magnitude of spin resistances. 
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